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Background and Purpose Variant alleles of CYP2C9 and VK0RC1 account for differences 
in anticoagulation response. We sought to establish a warfarin dosing formula for individu- 
alized target International Normalization Ratio of Prothrombin Times (INRs) using data from 
single nucleotide polymorphisms (SNPs) in VK0RC1 and CYP2C9 in Korean patients. 
Methods Ischemic stroke patients displaying stable target INR for at least 3 months before 
enrollment were analyzed. Warfarin and vitamin K levels were measured to adjust for con- 
founders. Phenotypes were defined using the 'warfarin response index' (WRI) defined as INR 
divided by the daily maintenance warfarin dose. We tested SNPs in CYP2C9 (3 sites: 430C>T 
(rs1799853), 1075A>C (rs1057910), 1076T>C) and VK0RC1 (14 sites: 381C>T, 861C>A 
(rs1 7880887), 2653G>C, 3673A>G, 5496G>T, 5808T>G (r1 78821 54), 6009C>T, 6484T>C 
(rs9934438), 6853C>T (rs1 7886369), 7566T>C, 8767G>C, 8814T>C, 9041G>A (rs1 7880624), 
and 9071 G>T) using a standard sequencing method. Multivariate linear regression analysis 
was applied to establish the formula for warfarin dosage. 

Results All 204 patients had excellent drug compliance. The mean INR was 2.22 (+0.56) 
and mean daily maintenance dose of warfarin was 3.92 mg (+1.54). Patients with low WRI 
were younger (P< 0.001) with high body mass index (P= 0.003), high prevalence of wild- 
type CYP2C9 polymorphism (1075A>C, P<0.001), and six heterozygote SNPs in VR0RC1 
(P<0.001), which were tightly interlinked (381T>C, 3673G>A, 6484T>C, 6853C>G. 
7566C>T, 9041G>A) (r 2 =1). Based on these data, a warfarin dosing formula was estab- 
lished. 

Conclusions WRI is influenced by age, body mass index and SNPs in VK0RC1 and CYP2C9 
in Korean stroke patients. The obtained warfarin dosing formula may be clinically applicable. 
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Introduction 

Warfarin is the most widely used anticoagulant in stroke pa- 
tients. Oral anticoagulation is superior to antiplatelet therapy 
involving clopidogrel or aspirin in the secondary prevention of 



stroke in patients with cardiogenic embolism. 1,2 However, war- 
farin has a narrow therapeutic index, particularly for stroke pa- 
tients; a suboptimal International Normalization Ratio of Pro- 
thrombin Time (INR) level would increase stroke recurrence 
while excessive anticoagulation results in brain hemorrhage. 
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Considering the high mortality of warfarin-induced brain hem- 
orrhages, proper adjustment of anticoagulation intensity is criti- 
cal for the effective management of stroke patients. 

Age, gender, concurrent medications, and dietary vitamin K 
intake influence the INR level, 3 7 and algorithms using these pa- 
rameters have been developed. 6,8 Nevertheless, these factors are 
insufficient to explain the wide inter-individual differences in 
doses required for optimal INR, 9 ' 11 which maybe related to in- 
dividual genetic variability. 11 ' 14 More recent studies have dem- 
onstrated that variant alleles of Vitamin K expoxide reductase 
complex, subunit 1 (VKORCl) and Cytochrome p450, sub- 
family IIC, polypeptide 9 (CYP2C9), the two enzymes related 
to the metabolism of vitamin K and warfarin, 7 ' 15 account for ge- 
netic variabilities in the anticoagulation response. 16 19 The inclu- 
sion of VKORCl and CYP2C9 to a dosing algorithm may en- 
able better prediction of warfarin dosing than when demogra- 
phic data alone are used. 20 22 Several warfarin dosage formulae 
and dose determination algorithms that include genetic infor- 
mation have recently been reported. 22,23 

The use of genetic polymorphisms in regulating warfarin 
dosase has recently been studied in Asian patients including 
Japanese, 24 ' 26 Chinese 27,28 and Koreans. 29 ' 31 However, these pre- 
vious algorithms have focused on the generalized warfarin dos- 
age adjustment without considering tailored individual target 
INR. Accordingly, we examined the single nucleotide polymor- 
phisms (SNPs) in VKORC 1 and CYP2C9 in Korean stroke pa- 
tients, and analyzed their effects on warfarin dosing to optimize 
stable INR. We additionally attempted to formulate a warfarin 
dosage equation that may be practically used for stroke patients. 

Methods 

Subjects 

Between November 2009 and August 2010 at the Asan Med- 
ical Center neurology outpatient clinic, we prospectively en- 
rolled patients who were diagnosed as having ischemic stroke 
confirmed by brain imaging studies, aged > 18 years, and ad- 
ministered warfarin regularly under strict instructions for sec- 
ondary prevention of ischemic stroke. Although INR level 2-3 
has been considered 'optimal' for secondary prevention of 
stroke, 32 considering its variation secondary to diet and occa- 
sional higher dosages needed in certain patients, such as those 
with mechanical valve replacement, 32 we considered the INR 
level 1.7-3.5 'relatively stable', 33 and enrolled patients who dis- 
played this INR level on more than three consecutive measure- 
ments over at least 3 months before enrollment. We excluded 
patients who were l) considered incompliant or anticipated to 
alter anticoagulation therapy (i.e., due to scheduled surgery); 2) 
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severely dementic or aphasic so that reliable communication 
was not possible; 3) taking long-term medication that could 
potentially influence anticoagulation intensity, such as anti-epi- 
leptics, antibiotics, anti-inflammatory drugs, and alcohol; 4) di- 
agnosed with underlying diseases potentially affecting the level 
of coagulation, such as vitamin K deficiency, malabsorption 
syndrome, malignancy; and 5) those who refused to consent. 
Patients with non-detectable R and S warfarin levels on the sam- 
pling date were additionally excluded. The study was approved 
by Ethics Committee of the Asan Medical Center, and written 
informed consent was obtained from all the subjects. 

We reviewed medical records and interviewed individual pa- 
tients. Data obtained were reasons for prescribing anticoagula- 
tion (e.g., arrhythmia, valvular heart disease), recent myocardiac 
infarction, heart failure, recurrent strokes despite antiplatelet 
use or severe intracranial stenosis, daily warfarin dose prescrib- 
ed, and the warfarin medication schedule. Demographic fea- 
tures and risk factors for stroke were also recorded including 
age, gender, hypertension (defined as receiving medication for 
hypertension or blood pressure > 140/90 mmHg with repeated 
measurements), diabetes mellitus (defined as receiving medica- 
tion for diabetes mellitus, fasting blood sugar > 126 mg/dL or 
post-prandial 2 hours > 200 mg/ dL with two separate measure- 
ments), hypercholesterolemia (defined as receiving cholesterol- 
reducing agents, overnight fasting low density lipoprotein (LDL) 
level > 100 mg/dL or overnight fasting cholesterol level > 200 
mg/dL), and current cigarette smoking (defined as smoking 
within six months of enrollment). The body mass index (BMl) 
was measured in each patient. 

Blood sampling and determination of biological 
phenotype 

Since 2000, we have provided patients receiving anticoagula- 
tion therapy official education with an instruction program. At 
the time of enrollment, a well-trained nurse confirmed that the 
patient had received the instructions. Blood samples (10 cc) 
were obtained 12 hours after warfarin intake. We measured INR 
using standard methods. 34 To determine patient drug compli- 
ance, the plasma warfarin levels of patients (R- and S -warfarin) 
were measured using chiral-phase high-performance liquid 
chromatography-tandem mass spectrometry (LC -MS/MS; 
ThermoFinnigan, San Jose, CA, USA), as described elsewhere, 
with minor modifications. 35,36 While all patients followed our 
diet education program, serum vitamin K levels (phylloqui- 
none) were additionally measured to adjust for the possible con- 
founding effects of diet on INR, using LC-MS/MS with atmo- 
spheric pressure chemical ionization, as described elsewhere. 37 

We assumed three variables as biological phenotypes influ- 
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enced by genotype: INR, daily maintenance warfarin dose, and 
warfarin response index (WRl). WRI was defined as INR di- 
vided by the daily maintenance warfarin dose (mg). Since WRI 
seems to be a more reliable variable reflecting individual re- 
sponse than INR or warfarin dose (see below), the index was 
employed as the biological phenotype in subsequent analyses. 
Based on WRI data obtained on the date of blood sampling, pa- 
tients were divided into three groups: high (upper, 33.3%), in- 
termediate (33.4-66.5%), and low (lower 33.3%) using two 
cutoff values. Novel SNPs were examined in two (high vs. low) 
groups. 

Determination of genotype 

We searched known SNPs related with warfarin metabolism 
from published papers, 31,38 database (http://www.ncbi/nlm. 
nih.gov/). We also sequenced whole VKORC1 genes for con- 
firmation of the presence of unreported SNP or linkage between 
the high WRI and low WRI groups. We selected 26 SNPs in 
CYP2C9 (three sites; 430C > T (rsl799853), 1075A> C (rslO 
57910), 1076T>C), VKORC1 (14 sites; 381C>T, 861C>A 
(rsl7880887), 2653G > C, 3673A> G, 5496G > T, 5808T > G 
(rl7882154), 6009C > T, 6484T > C (rs9934438), 6853C > T 
(rsl7886369), 7566T>C, 8767G>C, 8814T>C, 9041G>A 
(rsl7880624), 9071G > T), gamma-glutamyl carboxylase (GG- 
CX, nine sites; eight SNPs and one micro satellite; 268A>T, 
412G > A, 6258A> G, 8015C > T, 8016G > A, 8240A> G, 8445C > 
T, 13875A> G, CAA repeat). 

Sequencing method (resequencing analysis for the detection 
of genetic variations) was used to determine genotype. For se- 
quencing analysis, genomic sequences were obtained from the 
GenBank (http://www.ncbi/nlm.nih.gov/) database and PCR 
primers were designed using the Primer3 software (http: //www- 
genome.wi.mit.edu/cgi-bin/primer/primer3_www.cgi) to am- 
plify coding regions and promoter regions (> 500 bp upstream 
from exon l) ofVKORCl and CYP2C9 genes (PCR primer 
information used in this study is available on request). PCR 
amplification reactions were performed in 20-^L total volume 
containing 20 ng genomic DNA, 0.15 of each primer, 100 
uM of each dNTP, 1 unit of AmpliTaq Gold™ (Applied Biosys- 
tems, Foster City, CA, USA) and reaction buffer with 1.5 mM 
MgCL PCR was performed at 94°C for 5 minutes, followed by 
35 cycles at 94°C for 45 seconds, 55°C for 30 seconds, 72°C for 
1 minute and final extension at 72°C for 10 minutes. The PCR 
products (l ^L) were directly used as template for sequencing 
using BigDye™ terminator v2.0 cycle sequencing kit (Applied 
Biosystems) on an ABI3730 DNA sequencer according to man- 
ufacturer s instructions. DNA polymorphisms were identified 
using the PolyPhred program (http://droog.gs.washington. 
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edu/PolyPhred.html) after sequence chromatograms were base- 
called with the Phred program and assembled with Phrap. In the 
VKORC1, we searched unknown SNPs related with clinical 
phenotype. We selected high and low WRI groups and analyzed 
whole VKORC1 sequence to try and identify unknown SNP 
between groups. Also, we confirmed genotype using the Taq- 
Man SNP genotyping assay (Applied Biosystems) according to 
the manufacturer's instructions. Two polymorphic sites (rs7294 
for VKORC1 gene and rsl057910 for CYP2C9) were selected 
and genotyped for reconfirmation. Fluorescence intensities 
were read on a 7900HT Sequence Detection System (Applied 
Biosystems), and genotypes were determined using SDS 2.1 
software (Applied Biosystems). 

Statistical analyses 

Genotyping data were tested for deviation from the Hardy- 
Weinberg equilibrium using the Chi-square test. We analyzed 
log-transformed data of variables for biological phenotypes, 
specifically, INR, warfarin dose, and WRI, which were skewed. 
Fisher's exact test was used to compare the allelic frequencies of 
VKORC1 and CYP2C9 variants among patients with different 
variables. We calculated the effective sample size as 155.35 us- 
ing calculation fomula 39 (power 90%, a = 0.05, correlation coef- 
ficients = 0.1, Z (3 = 1.282, Z a/2 = 1.96, Fisher's transformed Z = 
0.5*ln (1.1/0.9) =0.1003, sample size=([l.96+1.282]/0.1003) 2 
+3 = 155.35). 

Relationships between items such as age, body weight, BMI, 
and WRI and warfarin, respectively, were analyzed using Pear- 
son correlation while the relationship between WRI and warfa- 
rin dose with items such as gender, hypertension, diabetes, hy- 
percholesterolemia, and cigarette smoking were tested with t- 
test. Multivariate linear regression analysis was performed to 
determine the formula for warfarin dosage. All statistical analy- 
ses were performed using SPSS 14.0 (SPSS, Chicago, IL), with 
the level of statistical significance being set at P < 0.05. 

Results 

Our study included 204 patients with a mean age of 66 years 
(SD, ±11.2) ranging from 25 to 87 years, including 128 males 
(63%). The mean BMI was 24.5 kg/m 2 (SD, ±2.8; range, 13.7 
to 32.2 kg/m 2 ). The reasons for anticoagulation included ar- 
rhythmia (150 patients, 73.4%), recurrent strokes despite anti- 
platelet use (39 patients, 19.2%), valvular replacement (10 pa- 
tients, 5%), and others (five patients, 2.4%). All patients showed 
good drug compliance. The mean INR was 2.22 (SD, ±0.56; 
range, 1.12 to 3.74), mean daily maintenance dose of warfarin 
was 3.92 mg (SD, ± 1.54; range, 1.29 to 10.0), and the mean 
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Table 1. Demographic features and risk factors, and correlation with biological 
phenotype (warfarin dose and WRI) 



Demographic features and 
risk factors (N= 204) 


Value 
(Mean±SD) 


Relation with 
WRI; Pvalue 


Relation with war- 
farin dose; Pvalue 


Age (year)* 


65.98+11.18 


< 0.001 


< 0.001 


Male(%) 


128(62.7%) 


0.194 


0.238 


Weight (kg)* 


65.1 ±10.54 


< 0.001 


0.001 


BMI(kg/mT 


24.45 ±2.80 


0.003 


0.006 


Hypertension 


83 (40.6%) 


0.294 


0.382 


Diabetes mellitus 


32(15.7%) 


0.423 


0.682 


Hypercholesterolemia 


68(33.3%) 


0.782 


0.924 


Smoking (%) 


3(1.5%) 


0.176 


0.229 



*P values according to Pearson correlation test, other items according to t-test. 



WRI was 0.648 (SD, ± 0.292; range, 0.178 to 1.806). 

As shown in Table 1, WRI was closely related to old age (P < 
0.001), BMI (P = 0.003), and body weight (P< 0.001 ). Mean 
R- and S-warfarin levels were 0.188 (SD, ± 0.1 10) and 0.081 ^g 
/mL (SD, ± 0.059), respectively and they were significantly cor- 
related to WRI (P< 0.001, respectively), while the vitamin K 
level (mean, 0.830 g/mL; SD, +0.581) was not (P = 0.123). 

We analyzed the prevalence of all targeted genotypes (data 
not shown). Seven SNPs were related to a significant difference 
of WRI (6 SNPs located in VKORC1 (38 IT > C, 3673G > A, 
6484T > C [rs9934438], 6853C > G [rsl7886369]. 7566C > T, 
9041G>A [rsl7880624]) and one in CYP2C9 [1075A>C]). 
In terms of genotype prevalence, 24 patients (11.8%) were het- 
erozygotic in six SNPs of VKORC1 and 185 patients (90.7%) 
were heterozygotic in one SNP of CYP2C9. The six SNPs in 
VKORC1 were 100% interlinked (r 2 = 1 ). In view of the finding 
that 9041 G>A (rsl7880624) displayed the most consistent 
results with repeated genotype tests using the TaqMan SNP ge- 
notyping assay (Table 2), we selected this site for subsequent 
analysis. 

The high prevalence of heterozygotic VKORC1 9041 G > A 
(rsl7880624) was closely associated with low WRI (P < 0.001), 
and prevalence ofheterozygote CYP2C9 1075A> C (rsl057910) 
with high WRI (P < 0.00 1) (Figure 1, Table 2). Other SNPs did 
not show notable correlation with WRI or other biological phe- 
notype variables examined. Moreover, in subgroup analysis, we 
sequenced the VKORC1 gene in patients with the high (upper 
33.3%) and low (lower 33.3%) WRI, but failed to identify novel 
SNPs that were significantly associated with biological pheno- 
type variables. Multivariate linear regression analysis identified 
four variables that independently influenced the WRI value: 
age, BMI, and two SNPs in VKORC 1 and CYP2C9 (R 2 = 0.275, 
multivariate linear regression, Table 3). Using regression coeffi- 
cients obtained from the above results, we established the war- 
farin dosage formula, as shown in Figure 2A. We additionally 
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Table 2. Characteristics of different genotypes 





VK0RC1 9041G>A 


CYP2C9 1075A>C 


Genotype 


(rs 7294) 


(rs 1057910) 






AA AC 


Prevalence 


180(80.2%) 24(19.8%) 


185(90.7%) 19(9.3%) 




[P= 0.948*) 


[P= 0.998*) 


Vitamin K 


0.837±0.593 0.775±0.487 


0.821 ±0.580 0.901 ±0.596 


(ug/mL) 


[P= 0.621) 


(P= 0.542) 


R- warfarin 


0.650 ±0.361 0.992 ±0.643 


0.713+0.416 0.474±0.377 


(ug/mL) 


[P< 0.001) 


(P=0.017) 


S-warfarin 


0.280 ±0.204 0.492 ±0.320 


0.309 ±0.230 0.273 ±0.230 


(ug/mL) 


[P< 0.001) 


(P=0.512) 


INR 


2.22±0.572 2.24±0.467 


2.24± 0.559 2.05± 0.555 




[P= 0.790) 


(P=0.155) 


Warfarin dose 


3.65±1.27 5.94+ 1 .83 


4.06±1.52 2.54±0.925 


(mg) 


[P< 0.001) 


(P< 0.001) 


WRI 


0.680 ±0.292 0.410 ±0.1 40 


0.625±0.278 0.877 + 0.331 




[P< 0.001) 


(P< 0.001) 



*P value in the Hardy-Weinberg equilibrium, Chi-square test; Univariate analysis, 
such as Student's t-test and Mann-Whitney test, with Bonferroni correction for di- 
chotomized variables. 



performed the above analysis using warfarin dose instead of 
WRI as a phenotype variable (Figure 2B). The results obtained 
were similar. 

Discussion 

We investigated the factors influencing WRI (and warfarin 
dose) in a relatively large number of homogeneous Korean 
stroke patients. Young age, high BMI, heterozygote VKORC 1 
(9041G>A, rsl7880624), and/or homozygote CYP2C9 
(1075A>C, rsl057910) contributed to low WRI, indicating 
that patients with these characteristics have low anticoagulation 
responses to a given warfarin dose. To minimize the confound- 
ing effect of non-compliance, we selected patients who had 
been formally educated with the institutional protocol. Further- 
more, since diet preference may influence the INR level, 40 we 
concomitantly measured the serum vitamin K level to adjust 
potential dietary effects on WRI. The vitamin K level did not 
influence the anticoagulation response. In addition, R- and S- 
warfarin levels were highly correlated with WRI (P < 0.00 1), il- 
lustrating the patients' excellent drug compliance. 

In this study, our aim was to determine the tailored intensity 
of anticoagulation response in Koreans and to provide a practi- 
cal warfarin dosage equation using WRI. The information about 
SNPs agreed with a previous genetic association study in Asian 
population showing the importance of genetic variability and 
prevalence of selected SNPs (VKORC1 9041G>A and CY- 
P2C9 1075A>C). However, the linkage in VKORC1 SNPs has 
not been discovered in the previous study in Korean popula- 
tion. 31 In the GGCX gene, we could not find any significant SNP, 
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Figure 1. WRI and warfarin doses in different genotypes. (A) WRI (left) and warfarin doses (right) in VK0RC1 rs7294; homozygote VK0RC1 rs7294 is associated with 
significantly higher WRI (P< 0.001 ) and lower warfarin dose (P< 0.001 ) than heterozygote. (B) WRI (left) and warfarin doses (right) in CYP2C9 rs1 05791 0; homozygote 
CYP2C9 rs1 05791 0 is correlated with significantly lower WRI (P< 0.001 ) and higher warfarin dose (P< 0.001 ) than heterozygote. 



Table 3. Contribution of age, BMI and selective genetic variations to biological 
phenotype using multivariate linear regression analysis 





Standardized regression 
coefficient (P value) for 
the estimation of WRI 


Standardized regression 
coefficient (Pvalue) for 
estimation of warfarin dose 


Age 


0.266 (< 0.001) 


-0.285 (< 0.001) 


BMI 


-0.132 (0.032) 


0.121 (0.034) 


VK0RC1 (rs 7294) 


-0.361 (< 0.001) 


0.413 (< 0.001) 


CYP2C9 (rs 1057910) 


-0.191 (< 0.001) 


0.276 (< 0.001) 


R 2 


0.275 


0.248 



In (INR/warfarin dose) 
= -0.395 + 0.01 (Age) - 0.02 (BMI) - 0.475 (V) - 0.279 (C 2 ) 



In (warfarin dose) 
= 1.142 - 0.01 (Age) + 0.017 (BMI) + 0.499 (V) + 0.370 (C 2 ) 



Figure 2. Warfarin dosage formula. (A) Formula using WRI (INR/warfarin 
dose) as a dependant variable. (B) Formula using the warfarin dose as a de- 
pendant variable. V, VK0RC; C, CYP2C9. 



which has shown inconsistent results in previous Japanese stud- 
ies. 24,26 

Our data were also consistent with results from studies with a 
multi-ethnic population 14,38,41 reporting that patients with at least 
one heterozygote in several SNPs in VKORC1 and CYP2C9 
1075A> C require a high warfarin dose. While 9041G > A was 
not linked with other SNPs (at positions 381, 3673, 6484, 6853, 
and 7566) among European-American, 38 our result in Koreans 
showed strong correlation (r 2 = 1 ) among six SNPs in VKORC 1. 



Also, we could not find significant influence of CYP2C9 (1076- 
T > C or 430C > T) on warfarin dose although this position 
was associated with warfarin dose in the previous studies. 38,41 " 43 
These differences maybe caused by ethnic genetic variability. 

In contrast to previous studies employing warfarin dose only 
as a phenotype variable 44 46 we used WRI, a variable that pro- 
vides a comprehensive estimation of the warfarin dose in rela- 
tion to target INR. While our main results are similar to those 
obtained using warfarin dose only, using WRI as a phenotype 
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variable has several advantages. First, the predictability of the 
formula using WRI (R 2 = 0.275) was higher than that with war- 
farin dose (R 2 = 0.248). Second, WRI distinguishes inter-indi- 
vidual biological phenotypes more precisely. For instance, pa- 
tients with INR 1.5 receiving 10 mg warfarin should be biologi- 
cally different from those with INR 4 receiving the same dose. 
Third, the formula using WRI as a variable is more adaptive in 
daily clinical practice. Since this formula includes the INR level, 
the warfarin dose can be adjusted according to the clinical need 
to alter target INR, such as when new ischemic or hemorrhagic 
strokes develop. Therefore, our data are practically applicable to 
general Korean population to calculate stable warfarin dose. 

Although our data may not be directly applicable to other 
ethnicities with different genetic background, our results may 
encourage further studies to develop the warfarin dosing for- 
mula in different ethnicities. 
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